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We present relativistic hydrostatic equations for obtaining irrotational binary neutron stars in quasi 
equilibrium states in 3+1 formalism. Equations derived here are different from those previously 
given by Bonazzola, Gourgoulhon, and Marck, and have a simpler and more tractable form for 
computation in numerical relativity. We also present hydrostatic equations for computation of 
equilibrium irrotational binary stars in first post-Newtonian order. 
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I. INTRODUCTION 

Preparation of reliable theoretical models on late in- 
spiraling stage of binary neutron stars is one of the most 
important issues for gravitational wave astronomy. This 
is because they are one of promising sources for grav- 
itational wave detector such as LIGO §, VIRGO @, 
GEO600 § and TAMA j§. From their signals, we will 
get a wide variety of physical information on neutron 
stars such as their mass, spin, and so on if we have a the- 
oretical template of them g|. In particular, a signal from 
very late inspiraling stage just prior to merging may con- 
tain physically important information on neutron stars 
such as their radius j^], which will be utilized for deter- 
mining equation of state of neutron stars || . 

Binary neutron stars evolve due to radiation reaction 
of gravitational waves, so that they never settle down to 
equilibrium states. However, the emission time scale will 
be always longer than the orbital period outside their 
innermost stable circular orbit (ISCO), so that we may 
consider that they are in quasi equilibrium states in their 
inspiraling phase even near ISCO. Motivated by this idea, 
there have been several works in which sequence of equi- 
librium states of binary neutron stars are computed and 
the sequence is regarded as an evolutionary track; for ex- 
ample, we have obtained corotational equilibrium states 
in first post-Newtonian approximation Q ; Baumgarte et 
al. have obtained corotational equilibrium states in a rel- 
ativistic frame work using conformal flat approximation 
g. Up to now, however, all relativistic works have been 
done assuming corotational velocity field |j| . As pointed 
out previously JlCfl , corotation is not a adequate assump- 
tion for velocity field of realistic binary neutron stars, 
because effect of viscosity is negligible for evolution of 
neutron stars in binary and as a result, their velocity field 
are expected to be irrotational (or nearly irrotational). 

For computation of realistic quasi equilibrium states 
of coalescing binary neutron stars just prior to merg- 
ing, Bonazzola, Gourgoulhon, and Marck (BGM) |ll]] 
recently presented a relativistic formalism. In their for- 



mulation, they assume a helicoidal Killing vector , and 
then project relativistic hydrodynamic equations onto a 
hypersurface orthogonal to I 11 . After that, they impose 
their irrotational condition on the hypersurface and de- 
rive hydrostatic equations for the irrotational fluid. We 
think, however, that there were several inadequate treat- 
ments in their work. First one is their definition of ir- 
rotational condition, because their irrotational condition 
is nothing but a necessary condition for irrotation even 
in the case when we assume existence of £ M |12|. In gen- 
eral case, their condition is not identical with the irro- 
tational condition. Second, in numerical relativity, we 
usually solve equations such as Hamiltonian constraint, 
momentum constraint, and equations for gauge condi- 
tions, using spatial coordinates on the hypersurface, £j, 
which is perpendicular to unit normal n v . Due to this 
reason, they had to re-project their equations onto £ t . 
As a result, their equations for determining velocity field 
have complicated form. Finally, in their formalism, it is 
necessary to solve a complicated vector Poisson equation 
for relativistic cases, which should be unnecessary for ir- 
rotational fluid. Although we may get correct results us- 
ing their formalism, we had better obtain a simpler and 
more tractable formalism. The purpose in this paper is 
to present such one. 

In section II, we derive hydrostatic equations for irro- 
tational fluid from relativistic hydrodynamic equations. 
We use 3+1 formalism and project the hydrodynamic 
equations onto Ej. Then, we impose an irrotational con- 
dition on E t , which agrees with the relativistic irrota- 
tional condition [^2|. As a result of projection onto E t , 
we obtain hydrostatic equations on E t , and hence, they 
have suitable forms to be solved in numerical relativity. 
Also, in our formalism, we need to solve only one Pois- 
son type equation for a scalar field for determination of 
vector field. In section III, taking Newtonian limit, we 
show that well-known Newtonian hydrostatic equations 
are derived from the present formalism. In section IV, 
we give first post-Newtonian hydrostatic equations for 
irrotational fluid as well as gravitational potentials to be 
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solved. Section V is devoted to summary. Throughout 
this paper, c denotes speed of light, and we use units in 
which gravitational constant is unity. We use units c = 1 
in section II for convenience and recover c in sections III 
and IV. Latin and Greek indices denote three dimensional 
(3D) spatial components (1 — 3) and four dimensional 
(4D) components (0 — 3), respectively. As spatial coordi- 
nates, we use the Cartesian coordinates x k = (x 1 , x 2 , x 3 ). 



II. RELATIVISTIC FLUID EQUATIONS IN 3+1 
FORMALISM 

Since we use 3+1 formalism in general relativity, we 
write the line element as 

ds 2 = g^dx^dx" 

= {-a 2 + f3 k (3 k )dt 2 + 2f3idx l dt + ^jdx'dx 3 , (2.1) 

where g^, a, = 7ij/3- J and 7 y are 4D metric, the lapse 
function, shift vector, and 3D spatial metric respectively. 
Using the unit normal to 3D spatial hypersurface St, 



and 



rip = I -a, 0,0,0 ), (2.2) 



jij is written as 



(2.3) 



Hereafter, we use and Di as the covariant derivatives 
with respect to g^ and 7 y , respectively. 

We assume the energy momentum tensor of the perfect 
fluid as 



P 



1 



P 

P 



u^ l u v + Pg>- 



(2.4) 



where p, e, P, and m m denote the rest mass density, 
specific internal energy, pressure, and four velocity, re- 
spectively. We assume polytropic equation of state P = 
(r — l)pe, where T = 1 + 1/n and n is the polytropic 
index. From adiabatic condition, we also get P = Kp T , 
where if is a constant. For the following, we define h as 



h= l 
= 1 



P 

■£ + — 
P 

r-r 



= i+/^. 

p 



(2.5) 



From the conservation equation for the energy momen- 
tum tensor, 



V„Tt = 0, 



we get the hydrodynamic equation as 
u^^Uu + V v h = 0, 



(2.6) 



(2.7) 



where u v — hu Ul and we use the conservation equation 
for rest mass density as 



V^pO =0. 



(2.8) 



To rewrite the hydrodynamic equation, we decompose 
it" as 



it" = u°(£> J - + V>*), 



(2.9) 



and assume that (1) is a timelike vector of its compo- 
nent (1, £ l ), and (2) U" is a spatial vector, U"n M — 0, i.e., 
V" = (0, V 1 ). By using £" and V", we get the following 
relations; 



7i 
7i 



Lf.U v - UyS v 
LlU v + /lV„ f-jr 



7 /£^ + ^A(^o) + {3) u k D t V k 

+na« <r 7<"^ i Vvn M , (2.10) 



and 



(3),~,. 



n <7 u <r 7 i , VV /l n„ ) (2.11) 



where £^ denotes the Lie derivative with respect to V 1 
and ( 3 )ui is a spatial vector d efine d as "f. k Uk- Using these 
relations, projection of Eq. (2.7) onto the 3D hypersur- 
face S t becomes 

tfCtUv + V k D k ^u l + ^u k D { V k + hDA\ 



+Dih = 0. 
We can rewrite this equation as 

+ V k {D k ^u l - A (3) «fc) = 0. 



(2.12) 



(2.13) 



Besides the conservation equation of the energy mo- 
mentum tensor, we ha ve the conservation equation for 
rest mass density (2.8). We note that for the ca se o f 
barotropic equation of state such as P = Kp r , Eq. (2.8) 
is also derived from the conservation equation of the en- 
ergy momentum tensor. This implies that if we solve Eq. 
( |2.8|) , we do not have to take into account n" component 
of Eq. Using Eq. (§J), E q- (fpl) is written as 



a[C e (pu°) + pu a V^] + D^pavPV 1 ) = 0. (2.14) 

Now, we assume that is a Killing vector such that 
V M C + Vvtf,. = 0, Ceu u = and Ce(pu°) = 0, and we 
write its component as (1, — tlx 2 , fix 1 , 0), where f2 is iden- 
tified with the orbital angular velocity with respect to 
distant inertial observer. We note that fluid exists in- 
side the light cylinder \x k \ <S cfi -1 , and existence of the 
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Killing vector is assumed within it. We also note that 
defined here is identical with the helicoidal Killing vec- 
tor defined by BGM |ll). If the Killing vector exists, 
we can derive hydrostatic equations for two interesting 
cases. One is the corotational case where we simply set 
V =0. Then, we get a well-known result as []l3f 



= constant, 



(2.15) 



and continuity equation is trivially satisfied in this case. 

The other is the case where satisfies an "irrota- 
tional condition" defined as 



and hence 



Di 



0, 



(2.16) 



(2.17) 



wher e (j) is a scalar field. Then, the hydrodynamic equa- 
tion (2.13) is integrated to give 



+ ^u k V k = constant. 



Note that V and u are written as 

v k = -e k -p k + -L^dkp, 



,o 



h-^D^Di 



1/2 



(2.18) 

(2.19) 
(2.20) 



so that we can rewrite left-hand side of Eq. ( 2.18| ) as 



± + &u k V k = ha 2 u° - (£ k + k )^u k . 



(2.21) 



By substituting Eq. ( |2.19| ) into Eq. ( 2.14 ), we get a 
Poisson type equation for determining <f> as 

Diipah^D^) - D^pau ^ + /?)} = 0. (2.22) 

Hence, hydrodynamic equations which should be solved 
for determination of equilibrium states reduce to only two 
hydrostatic equations ( 2.1 8| ) and ( 2.22| ). We do not have 
to solve any equations for vector potentials which were 
introduced in the formalism of BGM JO]. 

We note that definition of irrotation in the 4D covari- 
ant form should be |Ol 



= h^iy^Uy - V„%) = 0, 



(2.23) 



where PJ 1 = gj 1 + m ct u m , and we use Eq. (2.7) to rewrite 
the first line into the second line. When lo^ v is vanishing 
initially for a fluid element, it remains zero along the 
trajectory of the fluid element for the perfect fluid Jl2| . 
Hence, oj^ v — is just the irrotational condition. In 
our present irrotational condition (2.16), Eq. ( 2.23j ) is 
satisfied on the 3D hypersurface £t trivially. However, 



it is not trivial whether or not projection of Eq. ( 2.23| ) 
to n^^j v k component is satisfied. (Projection to n^n v 
component is trivially satisfie d.) W e here show that it is 
really guaranteed due to Eq. ( 2.16 ). By operating n^Y'k 
to V^u^ - V^Mp, we get 

n'V fe (V (tl u I , - V„f2 M ) 
= ~f u k C n u„ - 'f k (u M Vun» + n^V**^) 

= j v k C n u v + D k (hau°) ee W k , (2.24) 



where we use n^u. 



-hau°. From a straightforward 



calculation, we can rewrite 7 iy fc £„u„ as 



u r - 1 
7 k L n u v = - 

a 



Y k £iu„+hau Q D k a - (f3 j +P)d^u k 



(2.25) 



where d k denotes partial derivative on £ t . Hence, 



W k = - 

a 



«< 3 >fi* 



- (3) Ujd k (l3 j + £ j ) + d k (ha 2 u°) 



(2.26) 



Using the hydrodynamic equation ( 2.13| ) and an identity 
(2.21), we obtain 



W k 



-{V' + (3 j + ti){-d^u k + d k ^uj) 



(V j +p j +£ j )W kj . 



(2.27) 



Eq. (|2.27|) implies that W k = if Eq. fl2.16p is satisfied. 
Note that to derive Eq. fl2.27| ) we have not assumed the 
fact that £ M is a Killing vector. Therefore, Eq. ( 2.16 ) is 
the necessary and sufficient condition for the i rrotational 
condition in general case. Note that Eq. ( [2.18] ) itself does 
not mean irrotation in general. Even for the case when 
a Killing vector exists, it is nothing but a necessary 
condition for irrotation. 



III. NEWTONIAN LIMIT 

In the Newtonian limit, metric variables can be ex- 
panded as 



o = 1 - 5 + °( c ~ 4 )> 
f3 k = 0{c-% 

lij = $ij + 0(c~ 2 ), 



(3.1) 

(3.2) 
(3.3) 



where U denotes the Newtonian potential which satisfies 
AU = -471-p, (3.4) 
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and A is the flat Laplacian. By using v % = u z /u°, com- 
ponents of u M which we need here is also expanded as 



+ vA\ v2+u ) +0 ^ 



u l = Ui = h 0(c~ 3 ), 

c 



(3.5) 
(3.6) 



where v 2 = E^V. Note also ^ = (M7 C ) and V = 
(0,V' l /c). For the corotational case (V = ), we get 
Newtonian limit of left-hand side of Eq. ( 2.15 ) as 



full rcla. 



1 

1 + - 



v 2 f dP 



(3.7) 



Since V k — 0, v k is equal to £ k , and v 2 = R 2 Q 2 where 



B 2 



„1\2 , /„2^2 



(x ) . Substituting this relation of v into 



Eq. (3.7), we get a well-known result 



R 2 n 2 r dp 

U + — = constant. 

2 J p 



(3.8) 



For the irrotational case, Newtonian limit of the left- 



hand side of Eq. (2.1 q) becomes 



(3} u k V k 

2 



full rcla 



V 



-U 



f ^ + Y^v k (-e k + v k ) 

J P L. 



(3.9) 



In Newtonian order, v k = d k (j>-^, where <j> is expanded as 
<^>n/c+ 0(c~ 3 ). So that we get 

i ^(<9/c<M 2 - U+ j — £ k d k (t>N = constant, 

fe J P k 

(3.10) 



Eq. ( |3.10D agrees with that of BGM 

From continuity equations in Newtonian order, we ob- 
tain equations for 0n as 



(3.11) 



Eq. ( 3.11 ) is solved under boundary condition, 

^2(d k fa-t k )d k p = 0, (3.12) 



at stellar surface. 



IV. FIRST POST-NEWTONIAN EQUATIONS 

In this section, we derive hydrostatic equations in first 
post-Newtonian order. The equations for the corota- 
tional case agree with those shown in previous papers 



0, so that we here derive equations only for the ir- 
rotational case. In first post-Newtonian approximation, 
metric in the standard post-Newtonian gauge can be ex- 
panded as pa 



1 u 1 
a = 1 - — + — 



+ 0{c 



1 



c° 

lij = fiij 



k + O(c- 5 ), 



1 + %U 



+ 0(c~% 



where X and (3 k are obtained from 
AX = AnpUu + 2 ^(<9 fe N ) 2 + s 

^ k 

k = -\ph + \ (dk.x + ^^Pj 



ip 
p 



2 2 
and 

APfe = -47rp9jfe0 N , 
Ax = 47rp^(9fc^ N )a 



(4.1) 
(4.2) 
(4.3) 

(4.4) 
(4.5) 

(4.6) 
(4.7) 



Note that to derive these Poisson equations, we use a 
relation in Newtonian order, v k = d k (f>^. 
Using a post-Newtonian relation, 



au v = 1 



2c 2 



0n) 



1 



i(V(d fc N ) ; 



+ d k (t) N d k (j)p N - (rj + U) ^(a fc ^ N ) 5 

fe fe 
+ 0(c- 6 ), 



(4.8) 



where we expand <fi as <fo$/c + 0pn/c 3 + 0(c~ 5 ) an d ?/ = 
e + -P/p, first post-Newtonian expansion of Eq. (2. IS) 
becomes 



^u k V k 



full rcla. 



v-u + lY,( dk ^ 2 -J2 ekdk 

k fe 

-VU + \u 2 + X - ifo + 3C7) 5> fc N ) 2 

fe 



= constant. 



(4.9) 



First post-Newtonian expansion of continuity equation is 
also derived as 
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(4.10) 



where 

Ai = -t 



-rjdi<j)N - f3 t + di<j)PN 



(4.11) 



If is obtained from Eq. ( 3.11 ), Eq. ( 4.10 ) is regarded 
as equation for </>pn and solved under boundary condition 



A ^p = 



(4.12) 



at stellar surf ace. 

Eqs. (0|) and fl4.10| ) with Poisson equtions ([i~4"|), 
(4.6), (4.7) and (3.4) are basic equations for computation 



of irrotational equilibrium states in first post-Newtonian 
order. 



V. SUMMARY 

In this paper, we have derived relativistic hydrostatic 
equations for obtaining irrotational (quasi) equilibrium 
configurations of binary neutron stars using 3+1 formal- 
ism. In order to derive the hydrostatic equations, we 
first projected hydrodynamic equations onto Et and then 
impose the irrotational condition to obtain the hydro- 
static equations. As a result, the hydrostatic equations 
obtained are simple and suitable for numerical relativ- 
ity, compared with a previous formalism |Tl] ]. Also, as a 
natural consequence, in our formalism there is no vector 
Poisson equation to be solved and only a scalar Pois- 
son equation is needed to be solved for determination of 
velocity field not only for Newtonian, but also for rela- 
tivistic cases. 

We also give hydrostatic equations as well as Poisson 
equations for gravitational potentials needed for com- 
putation of irrotational equilirium states in first post- 
Newtonian approximation. We think that as a first step 
toward fully relativistic study, we had better construct 
post-Newtonian configurations for firm investigation of 
relativistic effect on binary neutron stars. In reality, we 
have been able to obtain much information on relativistic 
effect in binary neutron stars from first post-Newtonian 
studies |l6| llffl. Up to now, however, our attention 
was only paid to corotational binary neutron stars. The 
present formalism makes it possible to extend previous 
studies to the irrotational case. As a first work, we plan 
to obtain incompressible, irrotational equilibrium states 
of binary stars as we carried out for the corotational case 
previously |l6fj . 
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